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Abstract 

Background: Anopheles lesteri is a primary vector of Plasmodium spp. in central China. A complete understanding 
of vector population structure and the processes responsible for the differentiation is important to the vector- 
based malaria control programmes and for identifying heterogeneity in disease transmission as a result of discrete 
vector populations. There is no adequate An. lesteri population genetic data available. 

Methods: Polymorphism of sequence variations in mitochondrial COII and Cytb genes were assessed to explore 
the level of genetic variability and differentiation among six populations of An. lesteri from China. 

Results: There were 30 (4.37%) and 21 (5.33%) polymorphic sites for mtDNA-COII and Cytb gene, respectively. 
Totally 31 COII and 30 Cytb haplotypes were obtained. The range of fsr values was from 0.101 to 0.655 by mtDNA- 
COII, and 0.029 to 0.231 by Cytb gene. The analysis of molecular variance (AMOVA) showed that the percentage of 
variation within populations (65.83%, 88.48%) was greater than that among populations (34.17%, 11.52%) using 
both genes. The Tajima's D and Fu's Fs values were all negative, except Tajima's D values of YN and HNB 
populations, which suggest a large number of low-frequency mutations in populations and the populations were 
in expansion proceeding. 

Conclusions: Levels of genetic variation within An. lesteri populations were higher than among them. While these 
results may suggest considerable levels of gene flow, other explanations, such as the effect of historical population 
perturbations can also be hypothesized. 



Background 

Anopheles lesteri, which belongs to the Hyrcanus group 
of the genus Anopheles is a primary vector of malaria in 
central China [1]. Genetically-based methods have been 
proposed for malaria vector control. These methods 
focus mainly in altering vectorial capacity through the 
genetic modification of natural vector populations by 
means of introducing refractoriness genes or by sterile 
insect technologies [2]. Knowledge of the genetic struc- 
ture of vector species is, therefore, an essential require- 
ment as it should contribute not only to predict the 
spread of genes of interest, such as insecticide resistance 
or refractory genes, but also to identify heterogeneities 
in disease transmission due to distinct vector 



* Correspondence: yajunm@yahoo.com.cn 

Department of Pathogen Biology, Second Military Medical University, 
Shanghai 200433, China 



populations [3]. A complete understanding of vector 
population structure and the processes responsible for 
the distribution of differentiation is important to vector- 
based malaria control programmes and for identifying 
heterogeneity in disease transmission as a result of dis- 
crete vector populations [4]. Susceptibility to Plasmo- 
dium infection, survival and reproductive rates, degree 
of anthropophily, and the epidemiology of malaria in the 
human host may all be affected by genetic variation in 
vector populations [5]. 

Anopheles lesteri is almost morphologically undistin- 
guishable from its sibling species because of lacking the 
objective and stable identification characters, so the 
taxonomic status on An. lesteri in China has revised 
many times. Xu and Feng [6] regarded the Chinese "An. 
lesteri" as a new subspecies An. lesteri anthropophagus 
because it was distinct from both An. lesteri lesteri from 
the Philippines and An. lesteri paraliae [7] from 
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Malaysia in bionomics as well as morphology. The sub- 
species was later elevated to a full species rank [8]. 
However, the second internal transcribed spacer (ITS2) 
of ribosomal DNA (rDNA) of An. anthropophagus in 
China was similar to that of An. lesteri from the Philip- 
pines, South Korea, Guam and Japan [9,10]. The mole- 
cular evidence strongly support that An. anthropophagus 
is the synonym of An. lesteri. 

Anopheles lesteri exhibits variation in ecology [11], 
morphology [12], chromosomes [12], and random 
amplified polymorphic DNA (RAPD) markers [13]. 
Furthermore, An. lesteri was not considered as malaria 
vector in Guam and Philippines, but had high transmis- 
sion capacity of malaria in central China [11,14], and a 
certain transmission capacity in South Korea and Japan 
[15,16]. Despite its significance in malaria transmission, 
only a few studies on population genetics have been 
conducted [13]. Many genes of mtDNA were used to 
analyse the genetic variation and population structure of 
the Anopheline mosquitoes, such as cytochrome subunit 
I (COI) [17-20], cytochrome subunit II (COII) [21,22], 
control region [23], NADH dehydrogenase subunit 4 
[24] and subunit 5 [4,25-29]. The present study aimed 
to estimate genetic variability and population structure 
and to infer the extent of gene flow among An. lesteri 
populations from China based on mtDNA-COII and 
cytochrome B (Cytb) genes sequences. 

Methods 

Mosquito collections and species identification 

Wild adult An. lesteri were collected from 2004 to 2007, 
by using indoor light traps and human landing catches at 
human living room and livestock corrals. The eight col- 
lection sites in China were located from 22°17'N to 39° 
58'N, and 103°29'E to 123°50'E (Table 1 Figure 1). The 
HNB and YN populations consisted of specimens pools 
from two or three sites in proximity to each other, as sta- 
ted in Table 1. The distances between sites were below 
50 km. There were total five field populations and a 
laboratory colony, with JS population in this study. 



Adult mosquitoes of An. hyrcanus group were identi- 
fied by morphology using the identification keys of Lu 
et al [14]. Specimens were kept individually in silica gel 
filled tubes at 4°C, until DNA extraction was performed 
according to Collins et al [30]. Anopheles lesteri species 
identification was done by a PCR assay based on rDNA- 
ITS2 markers previously described in Ma et al [31]. 

mtDNA-COII and Cytb genes amplification and 
sequencing 

Sequence variation was examined in the mtDNA-COII 
and the Cytb genes. The COII and Cytb regions were 
amplified in 50 [iL reaction mixtures containing 1 x 
reaction buffer (QIAGEN, Courtaboeuf, France), 0.1 
mM of each dNTP (Eurogentec, Angers, France), 1 unit 
of Taq DNA polymerase, 0.1 |iM each of the forward 
and reverse primers and 1.5 \iL genomic DNA. The 
COII gene was amplified using primers COIIF (5'- TCT 
AAT ATG GCA GAT TAG TGC A -3', forward) and 
COIIR (5'- ACT TGC TTT CAG TCA TCT AAT G -3', 
reverse), and the Cytb gene using primers CytbF 
(5'-GGA CAA ATA TCA TTT TGA GGA GCA ACA 
G-3', forward) and CytbR (5'- ATT ACT CCT CCT 
AGC TTA TTA GGA ATT G -3', reverse). The cycle 
conditions in PTC- 100 Peltier Thermal Cycler included 
an initial denaturation step at 94°C for 2 min, followed 
by 30 cycles at 94°C for 30 s, 50°C for 30 s and 72''C for 
30 s, with a final extension step at 72°C for 8 min. After 
electrophoresis, PCR products were purified and used 
for sequencing in both directions with the previous pri- 
mers, on an ABI 3730 automatic sequencer (Applied 
Biosystems). Sequences were inspected and corrected, 
where necessary, using SEQSCAPE software (Applied 
Biosystems). 

Data analyses 

Multiple sequence alignments for each gene were per- 
formed using MEGA 4.0 [32] and CLUSTAL x [33]. 
The sequences polymorphism was assessed with MEGA 
4.0. A haplotype networks and outgroup probability of 



Table 1 The collecting data of Anopheles lesteri mosquito populations in this study 



Population Code 


Collecting site 


Collecting date 


Sample size 


Latitude (N) 


Longitude (E) 


GD 


Zhuhai, Guangdong 


Oct. 2007 


22 


22°17' 


113°30' 


YN 


Yanjing, Yunnan 


June 2006 


9 


28°60' 


104°13' 




Junlian, Sichuan 


June 2006 


4 


28°10' 


1 04°34' 


SC 


Pujiang, Sichuan 


June 2006 


23 


30°14' 


103°29' 


HNB 


Guangshui & Suizhou, Hubei 


June 2007 


8 


31"4V_3i°52' 


n3°15'-113°47' 




Tongbai, Heinan 


June 2007 


5 


32°29' 


113°23' 


JS* 


Xuyi, Jiangsu 


June 1985 


17 


32°54' 


n8°34' 


LN 


Donggang, Liaoning 


June 2004 


28 


39°58' 


123°50' 



* The collecting site and date of laboratory colony were original information. The mosquitoes were kept at 26 ± 1°C and 65 ± 5% (RH), under a 12: 12 hr (light: 
dark) photoperiod. 
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the haplotypes were constructed based on statistical par- 
simony using TCS 1.21 [34]. 

The parameters 0„ equivalent to the average pairwise 
number of differences between sequences [35], 6^ equiva- 
lent to the number of segregating nucleotide sites per 
sequence [36], and haplotypes diversity (h) were esti- 
mated for con and Cytb polymorphism within popula- 
tions. The population genetic structure was analysed with 
5 field populations, and assessed by analyzing molecular 
variance with ARLEQUIN 3.11 [37]. The percentage of 
sequence divergence within and between populations was 
calculated based on Nei and Li [38], and pairwise Fst 
values for short-term genetic distance between popula- 
tions were estimated with the methods of Slatkin (1995) 
[39] and tested for significance by permutation. Mismatch 
distributions were calculated using ARLEQUIN 3.11, and 
the neutrality tests were evaluated by Tajima's D and Fu's 
Fs. Isolation by geographical distance was assessed by 
GENEPOP 4.0.10 [40] using Mantel test. 



Results 

Sequences characteristics of mtDNA-COII 

One hundred and sixteen An. lesteri mosquitoes were 
distinguished by PCR assay from China (Table 1). A 686 
bp con sequence was determined in 88 mosquitoes, 
and a Cytb fragment of 394 bp was obtained from 112 
mosquitoes. All segregating sites and the sequence var- 
iants (haplotypes) are shown in Figures 2 and 3. The 
summary statistics for both genes are given in Table 2. 
Across the whole dataset, there were 30 (4.37%) and 21 
(5.33%) polymorphic sites for COII and Cytb, respec- 
tively. This low number of variable sites resulted in low 
nucleotide diversity and low haplotype diversity across 
samples. The 65 of overall field populations was from 
0.581 ± 0.435SD to 4.285 ± 1.709SD for COII, and 
0.274 ± 0.274SD to 3.545 ± 1.655SD for Cytb; 6^ was 
from 0.477 ± 0.485SD to 2.598 ± 1.606SD for COII, 
0.091 ± 0.188SD to 2.231 ± 1.476SD for Cytb and h was 
from 0.005 ± 0.003SD to 0.000 ± O.OOOSD (Table 2). 
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Haplotype 



Variable position GenBanl< 



con. 


_hl 


con. 


_h3 


con. 


.h4 


con. 


_h5 


con. 


.h6 


con. 


_h7 


con. 


_h8 


con. 


.h9 


con. 


.hll 


con. 


_hl2 


con. 


_hl3 


con. 


.hl4 


con. 


_hl9 


con. 


.h20 


con. 


_h21 


con. 


_h22 


con. 


_h23 


con. 


.h24 


con. 


_h25 


con. 


_h26 


con. 


.h27 


con. 


_h28 


con. 


_h29 


con. 


_h30 


con. 


.h31 


con. 


_h32 


con. 


_h33 


con. 


_h34 


con. 


_h35 


con. 


_h36 


con. 


h37 



1112223333444455555556666 Accession 
455683391792378034901237881278 „^ 

124184908515462406624727399204 ^_ 

GTGATGTCGTACTTTGCGATTAGCAATTTG EU699056 

A . T EU699058 

T T EU699059 

G. . . . EU699060 

T EU699061 

A . EU699062 

T C T EU699063 

T C EU699064 

TA T EU699066 

T C EU699067 

C T EU699068 

T EU699069 

.A C. . EU699074 

C . . EU699075 

T A EU699076 

T T EU699077 

. . . T G EU699078 

C EU699079 

A EU699080 

A EU699081 

. . T C EU699082 

..T C G... EU699083 

G EU699084 

T EU699085 

C EU699086 

T EU699087 

C T EU699088 

....C..T....C....A...G..G EU699089 

T A..G.... EU699090 

T....C G EU699091 

....C..T....CC G..G EU734846 

Figure 2 Variable bases of the mtDNA-COII gene for the haplotypes of Anopheles lesteri population 



Among the 88 COII sequences, 31 haplotypes were Cytb_2 and Cytb_4 were shared, especially; Cytb_2 

found. Four haplotypes of COII_l, COII_5, COII_6 and occurred in all populations (Table 2). Haplotype net- 

COII_20 occurred in more than one population, the fre- works showed that An. lesteri haplotypes derived from a 

quency was 12.90% (4/31). Thirty of 112 Cytb haplo- single common ancestral COII haplotype and two ances- 

types were observed. Three haplotypes of Cytb_l, tral Cytb haplotypes (Figure 4). 
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Variable position 



Haplotype 



3 1 2 3 5 7 8 



GenBanli 
Accession 
no. 





Z 


4 


i 


A 

4 


o 
0 


D 


a 


n 
I 


c 
D 


n 
U 


u 


4 


i 


n 

y 


1 


i 


/? 
D 


1 


i 


1 


J 




Cytb_l 


C 


G 


c 


T 


A 


A 


G 


T 


T 


G 


C 


T 


G 


A 


G 


T 


A 


C 


T 


T 


A 


EU699092 


Cytb_2 














A 






















A 








EU699093 


Cytb_4 














A 






A 
















A 








EUd99095 


Cytb_5 
























C 




















EU699096 


Cytb_6 












c 


A 






A 
















A 








EU699097 


Cytb_7 














A 






























EU699098 


Cytb_8 










G 




A 






A 
















A 








EU699099 


Cytb_10 














A 


















C 




A 








EU699101 


Cytb_ll 














A 




















G 


A 


G 


G 




EUd99102 


Cytb_12 










G 


G 


A 






A 
















A 








EU699103 


Cytb_13 








C 






A 




C 


















A 








EU699104 


Cytb_14 




A 










A 


C 




















A 








EUd99105 


Cytb_20 














A 


* 




• 










T 






A 








T7TTCAA •% •% ^ 

EU699111 


Cytb_21 














A 


c 




















A 








EU699112 


Cytb_24 














A 


c 












G 








A 








EU699115 


Cytb_26 














A 


















C 




A 




G 


G 


171 If /m -1 1 n 

EU699117 


Cytb_27 














A 


















C 




A 






G 


EU699118 


Cytb_28 






T 








A 


' 




• 












• 




A 








EUd99H9 


Cytb_29 








C 






A 


c 




• 








G 








A 








EU699120 


Cytb_3{J 


T 












A 


c 




















A 








1£Ud99121 


Cytb_31 


T 












A 






















A 








T7T IF* n n 4 n o 

EU699122 


Cytb_32 








C 






A 






A 
















A 








EU699123 


Cytb_33 








C 


G 




A 






A 












c 


G 


A 








EU699124 


Cytb_34 














A 








T 














A 








EU699125 


Cytb_35 








C 


G 




A 


















c 




A 








EU699126 


Cytb_36 








C 






A 






















A 








EU699127 


Cytb_37 














A 










C 












A 








EU699128 


Cytb_38 








C 


G 




A 






A 












c 




A 








EU699129 


Cytb_39 














A 












A 










A 








EU699130 


Cytb_4{) 








C 


G 




A 






A 


T 










c 




A 








EU699131 



Figure 3 Variable bases of the mtDNA-Cytb gene for the haplotypes of Anopheles lesteri population 



Population genetic structure of An. lesteri population 

The genetic structure was analysed with GD, LN, YN, 
SC and HNB populations. The range values of paiwise 
FsT was from 0.101 (GD/LN) to 0.655 (GD/SC) with 
mtDNA-COII, and 0.029 (HNB/LN) to 0.231 (YN/SC) 



with Cytb (Table 3). A Mantel test was carried out, and 
the correlation coefficient for the Fst with geographical 
distance was 0.271 by COII (P > 0.803) and 0.089 by 
Cytb {P > 0.400), which was not significance based on 
1,000 permutations. 
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Table 2 Data summary for populations, haplotypes and nucleotide diversity of Anopheles lesteri 



Pnni il;itinn 


Gene 


N 










5 


h 1+ 


SD) 

Jul 


e, (± SD) 


"n V— Jul 


VM 


Lull 


Q 
O 


1 ^/Tl Al~)\ 

\(Z), 5U), 4(/}, b[\} 








A 

4 


U.UUz ± 


U.UUz 


1 + 

\ .DHD JZ 


u.yo 1 


1 ^ -1- 1 TOT 

.043 ± 1 ./// 




LytD 


1 J 


1 'nc\ A('3\ c/i'i (^^^\ 7f'\^ 
\(Z), Z(b), 41 JJ, b[\), /[\) 








r 
J 


U.UU4 ± 


0.UU3 


1 f^l 1 + 
1 .O 1 1 IL 




.0 1 J ± 1.14/ 


MNd 


rr\\\ 
Lull 


Q 
O 


iM\ Q/^\ o/^\ ^^!'\^ I'i/'W i/in\ 








n 
o 


U.UU4 ± 


O.UUz 


Z.D \ 4 ± 


1 .dUo 


z.4C)4 ± 1 .oyz 




Cytb 


1 2 


1 n/i\ iA/i\ 1 

1 (3), 2(2), 4(1 ), 8(1 ), 1 0(1 ), 11(1), 1 2(1 ), 1 3(1 ), 1 


4(1) 






1 3 


0.006 ± 


0.004 


3.545 ± 


1.655 


l.li \ ± 1 .4/0 


JS 


coil 


12 


19(2), 20(10) 








1 


0.000 ± 


0.002 


0.331 ± 


0.331 


0.303 ± 0.379 




Cytb 


17 


2(17) 








0 


0.000 ± 


0.000 


0.000 ± 


0.000 


0.000 ± 0.000 


sc 


coil 


18 


6(3), 21(14), 22(1) 








2 


0.001 ± 


0.001 


0.581 ± 


0.435 


0.477 ± 0.485 




Cytb 


22 


2(21), 4(1) 








1 


0.000 ± 


0.000 


0.274 ± 


0.274 


0.091 ± 0.188 


GD 


coil 


18 


1(9), 23(1), 24(3), 25(1), 26(1), 27(1), 28(1), 29(1 








7 


0.002 ± 


0.001 


2.035 ± 


1.006 


1 .288 ± 0.946 




Cytb 


21 


2(6), 4(2), 20(1), 21(2), 24(2), 26(1), 27(3), 28(1) 


29(1 


, 30(1), 


31(1) 


10 


0.005 ± 


0.003 


2.780 ± 


1.239 


2.076 ± 1.350 


LN 


coil 


24 


1(4), 5(1), 7(1), 20(8), 30(1), 31(1), 32(3), 33(1), 


34(1) 


35(1), 36(1), 37(1) 


15 


0.005 ± 


0.002 


4.285 ± 


1.709 


2.598 ± 1 .606 




Cytb 


27 


2(14), 4(1), 32(1), 33(1), 34(1), 35(3), 36(1), 37(1 


), 38(2), 39(1 


),40(1) 


9 


0.005 ± 


0.003 


2.037 ± 


0.928 


1 .926 ± 1 .255 



h is haplotype diversity, 5 is the number of segregating sites and and Q„ are the estimates of nucleotide diversity. The bold haplotypes occur in more than one 
population and the number in parentheses indicates the frequency of the haplotype. 




Yang ef al. Malaria Journal 201 1, 10:216 
http://www.malariajoumal.eom/content/10/1/216 



Page 7 of 9 



Table 3 Pairwise genetic distance (Fsr) for populations of 
Anopheles lesteri 





YN 


HNB 


SC 


GD 


LN 


YN 


0.345/ 

0.114 


0.055 


0.231* 


0.128* 


0.153* 


HNB 


0.124 


0.324/ 

0.096 


0.125* 


0.055* 


0.029 


SC 


0.500* 


0.501* 


0.374/ 

0.159 


0.109* 


0.172* 


GD 


0.272* 


0.319* 


0.655* 


0.252/ 

0.099 


0.117* 


LN 


0.135* 


0.194* 


0.455* 


0.101* 


0.315/ 

0.103 



The pairwise values calculated by mtDNA-COII and Cytb gene are below and 
above the diagonal. The numbers along the diagonal are F^r values within 
population. The bold values are by COII gene. *P < 0.05. 



In the hierarchical AMOVE, both the 'among popula- 
tions' and 'within populations' variance components 
were considerable high, the latter was more contribution 
to total variances than the former (Table 4). The mean 
genetic divergence among populations was greater by 
COII (0.342) than Cytb (0.115). 

The simulated mismatch distribution among the 
mtDNA-COII and Cytb haplotypes was smooth and 
unimodal peak, which coincide with the population 
expansion model. Although, observed value appeared 
multimodal, the result of variance test indicated the 
degree of coincidence between them was not signifi- 
cance (P > 0.00 with COII, P > 0.15 with Cytb) [41]. 
The Tajima's D and Fu's Fs values were all negative, 
except Tajima's D values of YN and HNB populations 
(Table 5), which suggested a large number of low-fre- 
quency mutations in populations and the populations 
were in expansion proceeding. The strongly negative 
values for Fu's Fs suggested population growth and this 
is supported by the estimated values using COII gene 
from the rapid expansion model fitted in ARLEQUIN 
(t = 2\xt = 2.615, Bo = 0.00-0.39, Bj = 99 999, u = per 
sequence mutation rate, t = time since expansion, N = 
effective number of females). With a mutation rate of 
1 X 10'^ per site per generation [42], these values sug- 
gested a change in population size from a few thousand 
females to 10 females, in the range of 3970 years ago 
based on two generations of Anopheline mosquitoes in 
one month. 



Table 5 Values of neutrality test for Anopheles lesteri 



populations by mtDNA genes 







YN 


HNB 


SC 


GD 


LN 


Tajima's D 


COII 


0.283 


0.303 


-0.438 


-1.246* 


-1.402* 




Cytb 


0.009 


-1.498* 


-1.162 


-0.873 


-0.167 


Fu's Fs 


COII 


-8.139 




-28.504 


-26.580 


-26.580 




Cytb 


-18.044 


-15.265 




-27.052 


-27.112 



•P < 0.1. 



Discussion 

Sampling strategy and geographic coverage greatly influ- 
ence the analysis and interpretation of the data gener- 
ated from the samples. In China, An. lesteri was 
distributed in a range as the east of 100° E, and from 
19° N to 42° N [14]. In this study. An. lesteri mosquitoes 
were collected from most localities across its range. 
Although field An. lesteri specimen was difficult to col- 
lect due to usage of insecticide and environment 
changes, our sampling still covered geographic span of 
An. lesteri distribution. The LN was at the most north- 
ern limit, and GD was at the most southern limit of the 
distribution basically. 

In this study, both level of mtDNA- Oil and Cytb gene 
nucleotide diversity in field populations were greater than 
JS laboratory colony, such as all Cytb sequences in JS 
population were the same, which was similar to other 
gene on mitochondrial DNA, as COI {An. dirus, An. dar- 
lingi, An. stephensi) [17-20] and COII [An. jeyporiensis, 
An. minimus ) [21,22]. Thus, they are useful marker for 
exploring An. lesteri population genetic structure. 

The pairwise genetic distance using mtDNA-COII 
gene (0.101-0.655) was higher than Cytb (0.029-0.231). 
In theory, it was hard to prevent genetic divergence 
caused by genetic drift if the gene flow [Nm= (1- Fst)/^ 
Fst]) value was less than one [43]. The level of gene 
flow in these An. lesteri pairwise populations was below 
one, except YN/HNB, YN/LN, HNB/LN and LN/GD 
using mtDNA-COII gene, but all more than one except 
SC/YN using Cytb. The shallow population genetic 
structure was showed by Cytb gene. But the results by 
COII gene suggested that there was an apparent segre- 
gation from LN with the other populations, which is in 
agreement with the previous investigations with RAPD 
markers [13]. So, the level of An. lesteri population 



Table 4 AMOVA analysis of genetic variation in Anopheles lesteri populations by mtDNA genes 



Source of variation 


Degree of freedom 


Variance components 




Percentage of variation 




coil 


Cytb 


con 


Cytb 


COII 


Cytb 


Among populations 


4 


4 


0.433 


0.109 


34.17 


11.52 


Within populations 


71 


93 


0.835 


0.772 


65.83 


88.48 


Total 


75 


97 


1.268 


0.872 


100 


ICQ 



fsr = 0.342(COII), 0.115 (Cytb) 
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genetic divergence using mtDNA-COII gene should 
represent wild populations. 

The factors responsible for population genetic struc- 
ture should be analysed related with the climate, geo- 
graphy and the behaviour of mosquitoes. Yunnan is a 
highly complex region topographically due to its transi- 
tional position from tropical southern Himalayas to 
eastern Asia and from tropical Southeast Asia to sub- 
tropical China as well as at the junction of the India 
and Burmese plates, derived from Gondwanaland, and 
the Eurasian plate [44] . It is a noted centre of biodiver- 
sity [45-47]. It could have retained sufficiently mesic 
habitats for mosquitoes during the glaciations, when 
drier, more open habitats were spread widely [48]. If 
YN population of An. lesteri was the ancestor and the 
other region populations spread from Yunnan in the 
late stage of glaciations. The haplotype network sug- 
gested that An. lesteri migrated and spread from Yun- 
nan towards the North and the East China, and 
occurred colonization and expansion during migration 
proceeding. They were the same as the An. lesteri 
population patterns with An. dirus complex in South- 
east Asia by mtDNA-COI and microsatellite DNA 
[17,49], An. jeyporiensis in South China by mtDNA- 
COII [21]. If the migrating and expansion route was 
true, the An. lesteri samples in south of Yunnan should 
be increase to further investigation. An. lesteri is wide- 
spread in Palaearctic and Oriental region, and there is 
different climate, breeding habitation and blood prefer- 
ence, such as An. lesteri in southern and central China 
mainly is anthropophagic, but in Liaoning preferred 
animal's blood [11]. The above should be the key fac- 
tors of influencing population genetic structure of 
An. lesteri in China. 

Conclusion 

Levels of genetic variation within An. lesteri populations 
were higher than among them. There was an apparent 
segregation from Liaoning with the other populations 
using mtDNA-COII gene. The results of neutrality test 
suggested a large number of low-frequency mutations in 
populations and the populations were in expansion pro- 
ceeding. While these results may suggest considerable 
levels of gene flow, other explanations such as the effect 
of historical population perturbations can also be 
hypothesized. 
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